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We investigate the possibility that part of the dark matter is not made out of the 
usual cold dark matter (CDM) dustlike particles, but is under the form of a fluid of 
strings with barotropic factor w s = — 1/3 of cosmic origin. To this aim, we split the dark 
matter density parameter in two terms and investigate the dynamics of a spatially fiat 
universe filled with baryons, CDM, fluid of strings and dark energy, modeling this latter 
as a cosmological constant or a negative pressure fluid with a constant equation of state 
w < 0. To test the viability of the models and to constrain their parameters, we use the 
Type la Supernovae Hubble diagram and the data on the gas mass fraction in galaxy 
clusters. We also discuss the weak field limit of a model comprising a significant fraction 
of dark matter in the form of a fluid of strings and show that this mechanism makes it 
possible to reduce the need for the elusive and up to now undetected CDM. We finally 
find that a model comprising both a cosmological constant and a fluid of strings fits very 
well the data and eliminates the need of phantom dark energy thus representing a viable 
candidate to alleviate some of the problems plaguing the dark side of the universe. 
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1. Introduction 

Soon after the discovery of cosmic acceleration from the Hubble diagram of the high 
redshift Type la Supernovae (SNela) 1,2 , a strong debate arose in the scientific com- 
munity about the origin of this unexpected result. An impressive flow of theoretical 
proposals have appeared, while the observational results were constantly provid- 
ing more and more evidences substantiating the emergence of a new cosmological 
scenario. The anisotropy spectrum of the cosmic microwave background radiation 
(CMBR) 3 > 4 ' 5 , the matter power spectrum determined from the clustering proper- 
ties of the large scale distribution of galaxies 6 and the data on the Lya emitting 
regions 7 all provide indications that the universe have to be described as a spa- 
tially flat manifold where matter and its fluctuations are isotropically distributed 
and represent only about 30% of the overall content. In order to fill the gap and 
drive the acceleration, a dominant contribute from a homogeneously distributed 
negative pressure fluid has been invoked. Usually referred to as dark energy, the na- 
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ture and the nurture of this mysteryous component represent a new and fascinating 
conundrum for theoreticians. 

While the models proposed to explain this puzzle increase day by day, the most 
simple answer is still the old Einstein cosmological constant A. Although being 
the best fit to a wide set of different astrophysical observations 9,10 , it is neverthe- 
less plagued by two evident shortcomings, namely the cosmological constant problem 
and the coincidence problem. A possible way to overcome these problems invokes re- 
placing A with a scalar field (dubbed quintessence) evolving down a suitably chosen 
self interaction potential 11 . Although solving the cosmological constant problem, 
quintessence does not eliminate the coincidence one since too severe constraints on 
the potential seem to be needed thus leading to the fine tuning problem 12 . 

The ignorance of the fundamental physical properties of both dark energy and 
dark matter has motivated a completely different approach to the problem of cosmic 
acceleration relying on modification of the matter equation of state (EoS). Referred 
to as unified dark energy (UDE) models, these proposals resort to a single fluid 
with exotic EoS as the only candidate to both dark matter and dark energy thus 
automatically solving the coincidence problem. The EoS is then tuned such that 
the fluid behaves as dark matter at high energy density and quintessence (or A) at 
the low energy limit. Interesting examples are the Chaplygin gas 13 , the tachyonic 
field 14 and the Hobbit model 15 . 

It is worth noting that observations only tell us that the universe is accelerating, 
but they are not direct evidences for new fluids or modifications of the usual matter 
properties. It is indeed possible to consider cosmic acceleration as the first signal of 
the breakdown of the laws of physics as we know them. As a consequence, one has 
to to give off the standard Friedman equations in favour of a generalized version 
of them arising from some more fundamental theory. Interesting examples of this 
kind arc the Cardassian expansion 16 and the Dvali - Gabadadze - Porrati (DGP) 
gravity 17 both related to higher dimensional braneworld theories. In the same 
framework, one should also give off the Einsteinian general relativity and turn to 
fourth order theories of gravitation replacing the Ricci scalar curvature R in the 
gravity Lagrangian with a generic function f(R) that have been formulated both in 
the metric 18 ' 19 ' 20 and Palatini approach 21 < 22 < 23 > 24 providing a good fit to the data in 
both cases 25 > 26 . Actually, it is worth noting that it has been recently demonstrated 
that, under quite general conditions, it is possible to find a f(R) theory that predicts 
the same dynamics of a given quintessence model. 

Although resorting to modified gravity theories is an interesting and fascinating 
approach, it is worth exploring other possibilities in the framework of standard 
general relativity. Indeed, all the approaches we have described are mainly interested 
in solving the dark energy puzzle, while little is said about the dark matter problem. 
It is worth remembering that dark matter is usually invoked because of the need of 
a source of gravitational potential other than the visible matter. Considered from 
this point of view, it is worth wondering whether dark matter could be replaced 
by a different mechanism that is able to give the same global effect. Moreover, 
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such a mechanism must not alter the delicate balance between dark matter and 
dark energy that is needed to explain observations. Indeed, if we abruptly reduce 
the dark matter content of the universe without altering neither the dark energy 
term nor the background fundamental properties, we are not able to fit the available 
astrophysical data. Therefore, it is mandatory to test any proposed mechanism both 
at galactic and cosmological scales. 

In a series of interesting papers 28 , Letelier investigated the consequences of 
changing the properties of the right hand side of the Einstein equations adopting a 
fluid of strings as source term rather than the usual dust matter. Since such strings 
are not observed at the present time, it seems meaningful to extend the concept of 
dust clouds and perfect fluid referred to point particles to the case of strings. In 
particular, Letelier was able to find exact solutions for the case of a spherically sym- 
metric fluid of strings. It is worth noting that such strings could be of cosmological 
origin 29 and have thus to be included in the energy budget when investigating the 
dynamics of the universe. It is important to stress, however, that the strings we are 
referring to have finite lenght so that the results obtained for a network of cosmic 
strings of infinite length cannot be extended to the strings considered by Letelier. 
In this sense, the fluid of finite length strings we are considering represents a gen- 
eralization of the dust matter. While in this latter case, the matter particles are 
considered as pointlike, in the case of a fluid of strings 3, the elementary constituents 
are one dimensional objects with finite length. A fluid of strings has a profound 
impact at galactic scales. Indeed, assuming that the string transverse pressure was 
proportional to its energy density, Solcng 30 has demonstrated that the force law is 
altered thus offering the possibility of solving the problem of the flatness of spiral 
galaxy rotation curves 31 in a way similar to the MOND proposal 32 . 

Motivated by these considerations, we explore here the possibility that a part 
(if not all) of the dark matter may be replaced by a fluid of strings whose effective 
gravitational action may be considered as the source of the gravitational potential 
needed to flatten the rotation curves. To this aim, we consider cosmic strings as 
components of such fluid so that its e.o.s. may be simply parametrized by a constant 
barotropic factor w s = —1/3. Before discussing the impact at galaxy scales, it is 
preliminarily needed to investigate the effects at cosmological scales. We thus fit 
different cosmological models, both with and without such a component, to the 
available astrophysical data in order to test the viability of our proposal and explore 
if and how the constraints on the model parameters are affected by the presence of 
a fluid of strings. 

The paper is organized as follows. The models we discuss are described in Sect. 2, 
while the matching with observations is presented in Sect. 3 where we also compare 
the different models in terms of the information criteria parameters. Sect. 4 is de- 
voted to the weak energy limit of models comprising standard matter embedded in 
a fluid of strings and show how the corresponding modified gravitational potential 



a Hereafter, by fluid of strings we mean a fluid of finite length strings. 
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could help in reducing the need for CDM. A summary of the results and of their 
implications are presented in the concluding Sect. 5 



2. The models 

The key quantity entering most of the usual astrophysical tests is the Hubble param- 
eter H as a function of the redshift z. The position of the first peak of the CMBR 
anisotropy spectrum as measured by WMAP and balloon-borne experiments 3,4,5 
is a strong evidence of a spatially flat universe. Assuming therefore k — 0, the Fried- 
man equation for the expansion rate H = d/a (with a the scale factor normalized 
to 1 at today) reads : 

o 8irG 8irG v--^ /, % 

H = — ? T = — W 

where pi is the energy density of the i - th fluid and the sum is over the N cosmolog- 
ical fluids which make up the cosmic energy budget. If the fluids are not interacting, 
a conservation equation for each of them hold : 

Pi + ZH(l + w i )p i = Q (2) 

with the dot denoting derivative with respect to the cosmic time t and Wi = Pi/pi 
the barotropic facto of the i-th fluid. Assuming a cosntant Wi, Eq.(2) is easily 
integrated giving : 

p i (z) = Q i p crit (l + z)^ 1+w ^ (3) 

with Qi = pi(z — 0)/p cr it the present day density parameter of the z-th fluid and 
p cr it = 3i?o /8irG the present day critical density of the universe and, henceforth, 
we denote with a subscript all the quantities evaluated today. Inserting Eq.(3) 
into Eq.(l), we get : 



H{z) = H 0i 



N 



]T^(l + z)3(i+^) _ (4) 
\ i=i 

To fully assign the model, we have now to specify what are the ingredients of the 
cosmic pie and the values of their barotropic factors. According to the standard 
scenario, there are at least three components contributing to the energy budget, 
namely baryons, dark matter and dark energy. For the former two fluids, it is p = 0, 
while the dark energy is modelled as a negative pressure fluid with constant equation 
of state w < (as in quiessence models) with w = — 1 giving the usual A term. 
Motivated by the considerations discussed in the introduction, we add a fourth 
component 15 to our cosmological models. In order to see whether it is possible to 



b Note that we have not included radiation in the energy budget since its density parameter Sl ra( ; 



10 5 makes its contribute today indeed negligible. 
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reduce the dark matter content of the universe, we replace a fraction e of its energy 
contribute with a fluid of strings characterized by an equation of state 33 : 

w s = -1/3 (5) 

so that the energy density of the fluid of strings reads : 

p s (z) = n sPcrit {i + zf = eSl DM (l + zf (6) 

where hereafter the subscripts (b, DM, s, Q) denote quantities referred to baryons, 
dark matter, fluid of strings and dark energy respectively. The dimensionless Hubble 
parameter E{z) = H(z)/H thus finally writes: 

E 2 (z) = fi 6 (l + zf + (1 - e)JW(l + zf + 

+ en DM (l + z) 2 + n Q (l + zf 1+w ^ (7) 

where, because of the flatness condition, it is : 

Qq = 1 — Qf, — Vl DM . (8) 

Before investigating the consequences of introducing the fluids of strings, it is worth 
spending some more words on the philosophy underlying our model parametrization. 
Assuming that the dark matter is (mainly) made out of cold dark matter particles 
and denoting with the subscript CDM the related quantities, in Eq.(7), we have 
implicitly made the positions : 

&dm — ^cdm + 

(9) 

e = Q, s /VLt)m 

so that the parameter e gives an immediate feeling of what percentage of CDM 
may be given away without changing dramatically the dynamics of the universe, 
i.e. still obtaining a good fit to the available astrophysical data. In a sense, we are 
trying to reduce the need for dark matter replacing its contribute to the dynamics 
of the universe with a different kind of fluid having a different barotropic factor. 
Given our ignorance on the fundamental dark matter properties, there is no a priori 
reason against changing its equation of state. Moreover, it is also conceivable that 
the total dark matter turns to be made out of both CDM and the fluid of strings. 
As a consequence, the total matter density parameter is : 

fijVf = Slfe + flcDM + = + ^DM • (10) 

Eq.(7) refers to a cosmological model whose energy density is determined by four 
different components, namely baryons, dark matter, a fluid of strings and dark 
energy with constant w. Starting from this general case, we define four different 
models setting some of the parameters as follows : 



c Hereafter, we drop the subscript "Q" from us since this only refers to dark energy, while it has 
been fixed for both matter (baryons and CDM) and fluid of strings. 
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• ACDM: w = —l,e = 0. This is the popular concordance model that we 
consider as a testbed of our likelihood analysis. 

• ASDM: w = —1,0 < e < 1. Here, we still retain the cosmological constant 
as source of cosmic acceleration, but replace a fraction e of the dark matter 
term with the fluid of strings. 

• QCDM: w < — 1/3, e = 0. Also referred to as quiessence, this model repre- 
sents the easiest generalization of the successful ACDM scenario. Note that 
we do not impose a priori w > — 1 in order to explore the phantom models 
that seems to be favoured by the recent SNela data 2 . 

• QSDM: w < -1/3,0 < e < 1. This is similar to the QCDM considered 
above, but now we allow a fraction e of the dark matter to be replaced by 
the string fluid. A caveat is in order here. The fitting procedure does not 
choose a priori to decrease Sic dm while holding fixed Qq. As such, because 
of Eq.(8), it is possible that the search for the best fit ends in a region of 
the parameter space where CIq rather than Sic dm is reduced. 

As a general remark, we would like to note that although the four models above may 
formally be considered as a single one (since the former three are clearly particular 
cases of the latter one), they significantly differ in their underlying physics. As 
such, choosing among them is not only a matter of finding which one is in better 
agreement with the observations, but it is rather a sort of compromise between the 
capability of fitting the data and the physical justification of the model itself. 

3. Matching with observations 

Comparing model predictions with astrophysical observations is a mandatory test 
of the viability of the given model and also represents an efficient tool to constraint 
the characterizing parameters. We first describe the method we employ and then 
discuss the results of the fitting procedure. 

3.1. The method 

We fit the models described in the previous section to the SNela Hubble diagram 
and the data on the gas mass fraction in galaxy clusters. To take into account both 
datasets, we maximize the following likelihood function : 



£(p) oc exp 




(11) 



2 



with : 




(12) 
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Three terms enter the above \ 2 definition. The first one refers to the SNela Hubble 
diagram and is given by : 

NsNela r 

' / I i i — r I ; 

(13) 



XSNela — E 



where <7j is the error on the observed distance modulus ^i b s (zi) and the sum is over 
the number NsNeia of SNela observed. On the other hand, the theoretical distance 
modulus depends on the set of model parameters p and may be computed as : 



H{z) = 51ogL> L (z) + 25 
with Dl(z) the luminosity distance (in Mpc) given by: 

d( 



Dl{z) 



-(1 + *) 



H y 'Jo E((;p) 
The second term in Eq.(12) is defined as 34 ' 35 : 

fgas( z ii P) — fqas( Z i) 



v 2 
Agas 



E 



(14) 



(15) 



(16) 



with fgas(zi) the measured gas mass fraction in a galaxy cluster at redshift Zi and 
(Tgi the error. For a given cosmological model, f ga s(z,p) is given by 34 ' 35 : 

1.5 



fgas(z) 



M 



D s A CDM {z) 



DJ od (z) 



(17) 



where D A CDM and D™ od is the angular diameter distance for the SCDM and the 



model to be tested respectively. Da(z) may be evaluated as Da(z) = (1- 



2 D L (z). 



The constant b in Eq.(17) takes into account the gas lost because of different astro- 
physical processes. Following 36 , that have extensively analyzed the set of simula- 
tions in Rcf. 37 , we set b = 0.824. 

Finally, in Eq.(12), we have also included a prior on the age of the universe to 
that may be straightforwardly evaluated for a given set of model parameters as : 

dz 



to — tH 



I 

Jo 



(18) 



(l + z)E(z;p) 

where t H = l/H ^ 9.78ft, -1 Gyr is the Hubble time. 

With the definition (11) of the likelihood function, the best fit model parameters 
are those that maximize C(p). However, to constrain a given parameter pi, one 
resorts to the marginalized likelihood functions normalized at unity at maximum. 
The lcr confidence regions are determined by A% 2 — x 2 — Xo — 1> while the condition 
Ax 2 = 4 delimited the 2a confidence regions. Here, Xo i s the value of the x 2 f° r 
the best fit model. Projections of the likelihood function allow to show eventual 
correlations among the model parameters. In these two dimensional plots, the lcr 
and 2cr regions are formally defined by Ax 2 = 2.30 and 6.17 respectively so that 
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these contours are not necessarily equivalent to the same confidence level for the 
single parameter estimates. 

In order to reduce the space of parameters to explore, we fix the Hubble constant 
h (in units of 100 km s _1 Mpc -1 ) to the value determined by Daly & Djorgovski 
38 fitting the linear Hubble law to a large set of low (z < 0.01) redshift SNela: 

h = 0.664±0.008 . 

This value is in good agreement with h = 0.72±0.08 reported by the HST key 
project 39 based on a combined analysis of several local distance ladder methods. 
Since we arc not interested in constraining Hq, in the following analysis, we will set 
h = 0.664 neglecting the small uncertainty. Moreover, to a large extent, the effect 
of changing h on the results may be easily guessed and does not affect significantly 
the main results. 

The baryon density parameter f2{, is constrained by theoretical models of nu- 
cleosynthesis and by the observed abundance of light elements. Based on these 
considerations, Kirkman ct al. 40 have estimated : 

Q b h 2 = 0.0214±0.0020 . 

Combining this estimate with the value set above for the dimensionless Hubble 
constant h and neglecting the small error, we therefore fix = 0.049. 

3.2. Results 

We have applied the likelihood procedure described above using the SNela Gold 
dataset 2 and the catalog of relaxed galaxy clusters compiled in Ref. 36 . Moreover, 
we choose (tQ bs ,at) — (13.1,2.9) Gyr as obtained from globular clusters 41 and in 
agreement with estimates from nucleochronology 42 . 

Table 1. Summary of the results of the likelihood analysis of the models discussed in the text. The maximum 
likelihood value (bf) of each quantity is reported, while the 68% (95%) range is (bf — <5_ , 6/ + <5+) with <5— and 
<5+ the first (second) number reported as subscript and superscript respectively. The symbol (*) means that the 
parameter is held fixed. Note that, for the QSDM model, we may give only upper limits on e. 



Id 




e 




to (Gyr) 




ACDM 
ASDM 
QCDM 
QSDM 


n 970+0011 +0.016 
u._<u_ 002 _ .006 

n 970+0.007 +0.012 
U.Z(U_ ()004 _ .009 

n O7n+0 003 +0.008 
U.ZIU_ 0011 _ .016 

fl 97n+0-005 +0.009 
u - z,u -0.007 -0.015 


(*) 

n cn+0.15 +0.30 
U - OM -0.15 -0.35 

(*) 
(< 0.16 < 0.53) 


-!(*) 

-1 (*) 

i OO+0.07 +0.14 
i - z °-0.08 -0.15 

i 91 +0.10 +0.22 
1 " 21 -0.08 -0.17 


i o oy+0.02 +0.06 
±O.Z(_ o n _ .16 

ic t-i+0.59 +1.40 
lO.Ol_0.55 -1.09 

14 50+ ' 15 + - 29 
1 _.OU_ 12 _ .26 

14 60+ ' 38 + 1 ' 30 
1 *- ou -o .26 -0.47 


n R9+0.01 +0.02 
U - DZ -0.03 -0.04 

n Qfi+o- 19 +0- 26 

u - aD -0.11 -0.24 

n fin+0- 02 + 03 

U.OU _ 01 _ Q3 

n fi4+°- 06 +°- 21 

u - u -0.02 -0.05 



The results we get by applying the likelihood analysis presented above are re- 
sumed in Table 1 where we also give the estimated values of other physically inter- 
esting quantities, namely the age of the universe to and the transition redshift zt- 
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Fig. 1. Comparison among predicted and observed SNela Hubble diagram (upper panel) and 
fgas vs. z relation (lower panel) for the best fit ACDM model. 

This latter quantity is denned by the condition q{zr) = 0, being q = —aa/a 2 the 
deceleration parameter, and, for the general case of the QSDM model, it is given 
by: 

T ~[ n b + (i-e)n DM \ L - 

Since the uncertainties on the model parameters are not Gaussian distributed, a 
naive propagation of the errors is not possible. We thus estimate the 68% and 95% 
confidence ranges on the derived quantities by randomly generating 20000 points 
(Qdm,w,s) using the marginalized likelihood functions of each parameter (if not 
held fixed) and then deriving the likelihood function of the derived quantity. Al- 
though not statistically well motivated, this procedure gives a conservative estimate 
of the uncertainties which is enough for our aims. Let us now briefly discuss the 
results for each model. 



3.2.1. ACDM 

Not surprisingly, the ACDM model gives an almost perfect fit to the dataset con- 
sidered, as shown in Fig. 1. Having set from the beginning the Hubble constant, 
the only quantity to be determined is the dark matter density parameter Qdm 
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that turns out to be severely constrained by the likelihood test. Adding to ft dm 
the baryon contribution Q^ as set above, we get Qm = 0.319 as best fit value, 
while the 95% confidence range is (0.313,0.335). It is worth comparing our result 
with those previously obtained by other authors. For instance, using only the SNela 
Gold dataset, Riess et al. 2 (hereafter R04) have found Qm = ®-29-om- ^ n the 
framework of the concordance model, a combined analysis of the CMBR anisotropy 
spectrum measured by WMAP, the power spectrum of SDSS galaxies, the SNela 
Gold dataset, the dependence of the bias on luminosity and the Lya power spec- 
trum lead Seljak et al. 10 (hereafter S04) to the estimate Q M = 0.284±°;°™ (at 99% 
CL). Finally, fitting to the f gas data only with priors on both h and Qbh 2 , but not 
imposing the flatness condition ab initio, A04 estimates Qm = 0.245^q 037 : while 
including the CMBR data, they get Qm = 0-26to.o4- Overall, there is a very good 
agreement with our result. Nevertheless, it is worth noting that our best fit value is 
sistcmatically larger than that commonly quoted. This is partly due to having set 
Qb — 0.05 which is slightly larger than the fiducial value Qb ~ 0.04 often adopted. 
If we had set h = 0.72, the best fit value for Qm should be lowered by ~ 0.01 thus 
further reducing the difference with the standard result Qm — 0.3. 

It is worth stressing that the substantial agreement among our estimated Qm 
and the previous results obtained using a variety of methods makes us confident 
that the likelihood analysis we have performed is correct and is not affected by 
some systematic errors. It is thus meaningful to apply this method to the other 
models presented in Sect. 2. 

Although not directly constrained by the fitting procedure, it is nonetheless in- 
teresting to compare the derived quantities reported in Table 1 with other estimates 
in literature. First, we consider the age of the universe to whose best fit value turns 
out to be 13.27 Gyr. The 95% confidence range is within the prior set on this quan- 
tity, but, as we will see later, this is not a general result. Most of the more recent 
estimates of to have been obtained as a byproduct of fitting the ACDM model to 
a combination of different datasets and are thus rigorously model dependent. For 
instance, Tegmark et al. 9 (hereafter T04) give to = 13-24^ ; 41 Gyr, while Rebolo 
et al. 5 find t = 14. 4+^3 Gyr (all at la level). Both these estimates agree with 
to = 13.6±01.19 Gyr obtained by S04 that is the most comprehensive analysis. The 
accordance of our estimated to with these results is not surprising due to the fact 
that we are using the same model and have yet obtained a similar value for Qm- 

It is therefore more interesting to consider the constraints on the transition 
redshift zt- Since, for a flat ACDM model, zt only depends on Qm, constrain- 
ing zt is equivalent to constrain Qm- Unfortunately, it is difficult to derive this 
quantity directly from the data even if some attempts have been made. Fitting the 
phenomenological parametrization q{z) = qo + dq/dz\ z=0 z to the SNela Hubble dia- 
gram, R04 have found zt = 0.46±0.13 (at la level) which is in marginal agreement 
with our 95% confidence range. Since it is not clear what is the systematic error in- 
duced by the linear approximation of q(z), which only works over a limited redshift 
range, we do not consider a serious shortcoming of the model the la disagreement 
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Fig. 2. The same as in Fig. 1 but for the ASDM model. 

between the R04 estimate of zt and the one reported in Table 1. 
3.3. ASDM 

Let us now consider the results obtained for the ASDM model in which the cos- 
mic acceleration is still driven by the cosmological constant A, but the total dark 
matter content is made out by dust - like CDM particles and a fluid of strings. As 
clearly shown by Fig. 2, the model is able to fit very well both the SNela Hubble di- 
agram and the f gas data. Moreover, the estimated SJm (and thus CIm) is in perfect 
agreement with that obtained for the ACDM model and hence with all other re- 
sults discussed before. This could be qualitatively explained by noting that a fluid 
of strings is unable to drive cosmic acceleration even if it has negative pressure. 
Therefore, the amount of dark energy needed to accelerate the universe is the same 
as in the case of the ACDM model so that = 1 — CIm (and hence Hm) must be 
the same. 

The most striking result is, however, the constraint on e, i.e. the fraction of dark 
matter represented by the fluid of strings. The best fit value turns out to be 0.59 and, 
what is more important, the value e = is safely excluded at more than 2a level. 
Looking at Fig. 3, where the likelihood contours in the plane (flnM,£) are plotted, 
shows that e and SIdm are positively correlated (even if weakly) so that it is not 
possible to reduce e without decreasing the total dark matter content. We may thus 
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Fig. 3. 68%, 95% and 99% confidence contours in the (Qdm,e) plane for the ASDM model. 

safely conclude that is possible to fit cosmological data with a significant fraction of 
the dark matter content in the form of a fluid of strings rather than CDM particles. 
This conclusion is reflected in the estimated ilf, + flcDM that turns out to be much 
smaller than in the ACDM case, the best fit value being 0.154. This result may have 
interesting implications. Let us remember that a possible method to estimate 
consists in estimating the mass - to - light ratio M/L of clusters of galaxies and then 
integrating over the clusters luminosity function. Applying this method usually gives 
Qm — 0.16 (see, for instance, 43 and references therein) in striking disagreement 
with the results from tests probing cosmological scales (as SNela Hubble diagram 
and CMBR anisotropy spectrum). If we assume that the strings have a negligible 
mass d (which is a reasonable hypothesis 33 ), we may qualitatively conclude that 
the method outlined above should give an estimate of fit, + Clc dm rather than 
since it is unable to weigth the contribution of the fluid of strings. Actually, the 
high value of e found has profound implications also at a galactic scales as it will 
be discussed in much detail later. 

Although the picture of the today universe is similar to that of the ACDM model 
(i.e. the A term dominates over the matter one), the presence a non - negligible 
fluid of strings alters the dynamics of the universe introducing, for instance, a fluid 
of strings dominated period. As a consequence, the age of the universe (which is 
sensitive to the full evolutionary history) is significantly changed with the best fit 



d Note that this by no means imply that Q s is negligible since this is an energy rather than a mass 
density parameter. This could be best understood considering the case of radiation. Photons have 
zero rest mass, but nonetheless Q ra d does not vanish today and was dominant in the past. 
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value to = 15.51 Gyr that is outside the 95% range obtained for the ACDM model. 
Moreover, this values is also larger than our best fit prior t = 13.9 Gyr based on 
globular clusters and therefore model independent. However, to depends linearly 
on h~ x so that it is possible to decrease its value by simply increasing h without 
changing the other model parameters (and still obtaining a very good fit to the 
data). For instance, using h — 0.72 gives to = 14.3 Gyr in good agreement with 
the Rebolo et al. 5 estimate quoted above. On the other hand, to should be used 
to discriminate between the ACDM and ASDM model since they equally fit the 
same dataset, but predict significant different values for this quantity. Therefore, an 
accurate and model independent estimate of the age of the universe should make it 
possible to conclusively select one of the two models. 

The presence of the fluid of strings also affects the transition redshift which turns 
out to be much higher (0.96 vs. 0.62 for the best fit values) than in the ACDM case 
(and hence more in disagreement with the estimate of R04). This result could be 
qualitatively explained considering that, in order to have deceleration, the universe 
must be dominated by dust matter. Therefore, introducing a non - negligible fluid of 
strings component, delays the onset of dust matter domination and thus increases 
zt- If a reliable determination of this parameter were available in the future, we 
should obtain a further tool to confirm or reject the presence of the fluid of strings. 



3.3.1. QCDM 

Relaxing the hypothesis w = — 1 but still keeping it constant and assuming again 
that there is no fluid of strings, we get the QCDM model where there are two 
parameters to be constrained, namely the dark matter density parameter £Idm and 
the barotropic factor w. The best fit is obtained for (CIdm, w) = (0.270,-1.28) and 
it is shown in Fig. 4, while we refer to Table 1 for the constraints on the parameters. 
In particular, we note that the results on Qdm (and hence on Qm) are in perfect 
agreement with those obtained for the ACDM model so that we do not discuss 
anymore this parameter. 

It is more interesting to look at the constraints on the barotropic factor. First, 
we note that values of w < — 1 are clearly preferred, while w > — 1 is excluded at 
more than 95% level. In particular, the cosmological constant is ruled out by the 
likelihood analysis. Moreover, Fig. 5 shows that w and O^m are positively correlated 
so that increasing to is only possible by unrealistically increasing the matter con- 
tent. Surprising as it is, this result is however in agreement with previous analyses. 
Combining WMAP anisotropy spectrum with large scale structure clustering data 
and an old compilation of SNela, Spergel et al. 4 have found w — — 0.98±0.12 when 
dropping the prior w > — 1. Repeating the same analysis but using the Gold SNela 
sample, R04 have found w = — 1.02+ ; 19 . T04 added to the above dataset the power 
spectrum determined from the SDSS galaxy sample thus deriving to = — 0.721q; 3 2 7 . 
Finally, the comprehensive analysis of S04 gives w — — I.O8O+0.198 ( an results given 
at la). All these estimates agree among each other and our one, although we note 
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Fig. 4. The same as in Fig. 1 but for the QCDM model. 

that our best fit value is significantly smaller. Moreover, our result is the only one 
excluding the ACDM model at more than 99% level considering the data which 
have been used. 

Having pushed downward the confidence range for w affects the predicted age 
of the universe. The maximum likelihood value turns out to be to = 14.50 Gyr, 
while, at 95% level, to hes between 14.24 and 14.79 Gyr. These values turns out 
to be higher than those for the ACDM model, but lower (even if in agreement 
at the 95% level) than what is predicted in the ASDM case. As yet noted above, a 
comparison with the results obtained for the QCDM model by T03 and S04 suggests 
that including other kind of data (namely the CMBR anisotropy spectrum and the 
galaxy power spectrum) pushes upwards the constraints on w with values closer to 
the ACDM ones. A similar conclusion also holds for the age of the universe. For 
instance, Tegmark et al. give to — 13.53lo'g5 which is, however, in agreement with 
our estimate. Note also that our value could be reconciled with the T03 estimate 
by increasing the value of h from 0.664 to the value 0.71 used by these authors. 

As a final remark, we note that the transition redshift Zt is only weakly affected 
by relaxing the hypothesis w = —1 as it is witnessed by the good agreement between 
the estimates reported for the two cases. Indeed, zt mainly depends on the balance 
between and fig so that, being Qdm and the same in the two models, the 
resulting zt are naturally concordant. 
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Fig. 5. The same as in Fig. 3 but for the QCDM model. 

3.3.2. QSDM 

Finally, we discuss the results for the general QSDM model where we relax the 
hypotheses on both the barotropic factor w and the fraction e of dark matter made 
out of fluid of strings. Having now three parameters to constrain (SIdm, w,e), it is 
not surprising that the confidence ranges enlarge. While this is only a minor effect 
for what concerns SIdm, there is a significant weakening of the constraints on w and 
a dramatic impact on e on which we are only able to give upper limits. Nevertheless, 
the results are quite interesting. In particular, the best fit is obtained for a model 
that is very similar to the best fit QCDM one with the same content of dark matter 
(i.e. the same Qdm), almost the same w (-1.21 vs. -1.28) and no fluid of strings. 
Not surprisingly, the fitting to the data is perfect as shown in Fig. 6. 

Although quite similar to the QCDM model for the best fit values, the QSDM 
case is however different for what concerns the constraints on its parameters because 
of the presence of a third quantity to be determined. This is pictorially shown in 
Fig. 7 where we plot the projections of the likelihood function on the three orthog- 
onal planes in the parameter space. It is worth noting that £Idm is the parameter 
least affected by the addition of a third quantity so that the resulting constraints are 
only marginally weakened and the 68% and 95% confidence ranges nicely overlap 
with those obtained for all the models considered above. As an obvious consequence, 
our estimate for VIm (that is, Clb + flDM) is also in agreeement with the other results 
quoted above. Much care is needed when considering the constraints on the derived 
fib + flcDM- Since the best fit model has a vanishing e, the maximum likelihood 
value turns out to be almost the same as those for the ACDM and QCDM model 
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Fig. 6. Same as Fig. 1 but for the QSDM model. 




Fig. 7. The same as in Fig. 3 but for the QSDM model. In each panel, the value of the parameter 
not shown on the axes is set to its best fit value reported in Table 1. 

that also have no dark matter in fluid of strings. Nevertheless, the upper limits on e 
leads to strongly asymmetric constraints on the standard matter density parameter 
ilb + Qcdm so that, at the 95% level, also values as low as 0.166 are possible thus 
being more in line with the prediction of the ASDM model. 

As a general rule, the QSDM model is indeed very similar to the QCDM one, 
having almost the same matter content (and hence the same dark energy content) 
and constraints on the barotropic factor w which agree very well with those of the 
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Table 2. Number of parameters k, 
AIC A and BIC B for the models dis- 
cussed in the text. 



Id 


k 


A 


B 


ACDM 


1 


225.687 


228.896 


ASDM 


2 


215.830 


222.249 


QCDM 


2 


210.933 


217.352 


QSDM 


3 


213.456 


223.084 



QCDM model, but are less stringent. As a result, we get also concordant estimates 
of both the age of the universe and the transition redshift so that we do not comment 
again on these quantities. 



3.4. Comparing the models 

The results of the likelihood analysis discussed above have shown that all the models 
we have considered are able to fit well the same dataset. The next natural step is 
wondering which is the better one. To answer this question is not an easy task. 
Combining different datasets requires the introduction of the pscudo-x 2 defined 
in Eq.(12) that is not the same as the x 2 commonly used in statistical analysis. 
Moreover, the models have a different number of parameters to be constrained. As 
such, it makes no sense comparing the models on the basis of the \ 2 best fit value. To 
overcome this difficulty, Liddle 44 have proposed to resort to the information criteria 
that are widely used in other branches of science (such as medical pathologies), but 
poorly known in astrophysics. In particular, Liddle have proposed to use the Akaike 
information criterion (AIC) defined as 45 : 

A = -21n£ + 2fc (22) 

and Bayesian information criterion (BIC) 46 : 

B = -2\n£ + k\nN (23) 

with C the likelihood evaluated for the best fit parameters, k the number of model 
parameters and N the total number of points in the dataset used. Since both AIC 
and BIC explicitly takes into account the number of parameters, it is meaningful 
to compare the models on the basis of the values of these quantities. The lower is 
A or B, the better the model is. Unfortunately, as discussed by Liddle, it is not 
an easy task to decide what is the better information criterion to be adopted so 
that we conservatively report in Table 2 the values of both A and B for the best fit 
model parameters. The models are quite similar in their AIC and BIC values, as it is 
expected since they equally fit the same dataset. Nevertheless, it is possible; to rank 
them according to A or B. Independently on what criterion is used, the QCDM 
model turns out to be the preferred one followed by the QSDM or the ASDM 
depending on which information criterion (AIC or BIC respectively) is adopted. 
This simple comparison leads to two quite interesting conclusions. 
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First, there is a clear evidence favouring models others than the concordance 
ACDM. In particular, the AIC suggests that phantom like models have to be pre- 
ferred so that a violation of the strong energy condition is unavoidable to explain 
the cosmic acceleration in the framework of constant w dark energy models. 

A second and perhaps more interesting result is that the ASDM model is pre- 
ferred over the ACDM model and only marginally disfavoured with respect to the 
QCDM one. Therefore, introducing a fluid of strings as a component of the dark 
matter term makes it possible not only to better fit the data, but it could also be 
a viable alternative to phantom like models. Surprising as it is, this result is quite 
encouraging and motivates further study. 

4. The weak field limit 

The results discussed above has shown that introducing a fluid of strings in the 
dark matter budget does modify the dynamics of the universe, but gives rise to 
a model that is still in agreement with the SNela Hubble diagram and the data 
on the gas mass fraction in galaxy clusters and predict an age of the universe 
which is not unreasonable. Moreover, both the Akaike and Bayesian information 
criteria quantitatively indicate that the ASDM model has to be preferred over the 
concordance ACDM. 

It is worth stressing that, for the ASDM model, the fraction of dark matter 
represented by the fluid of strings is nearly dominant so that the density parameter 
of the standard (baryons + CDM) is significantly smaller with a best fit value as 
low as 0.154. It is easy to understand that lowering SIcdm have a profound impact 
at galactic scales. Indeed, since we have assumed that the strings have a very small 
(if not vanishing) mass, a small flcDM automatically implies less massive dark 
matter haloes. It is thus worth wondering whether such light haloes may still fit 
the rotation curves of spiral galaxies. Naively, one should think that the answer is 
negative since values of v c (r) larger than those predicted on the basis of luminous 
matter only naturally invoke massive haloes. But this is only true in a Newtonian 
gravitational potential. This is no more the case in a ASDM model. Considering the 
weak field limit, Soleng has shown that the gravitational potential for a pointlikc 
mass m embedded into a halo of fluid of strings is given by° 30 : 



where c is the speed of light, £i ; 2 are two integration constants and a = —p/p. 
Since, for our model, it is a = 3, imposing the condition that $(r) reduces to the 
usual Newtonian potential in the case r << £1,^2, we get £1 = 2Gm/c 2 so that, 




(24) 



°It is importnat to stress that we are considering finite length strings rather than a network of 
cosmic strings with infinite length. Indeed, in this second case, it is still not clear how to compute 
the gravitational potential in the weak field limit. 
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with simple algebra, we can rewrite Eq.(24) as: 



$(r) = - 



Gm 



1/3" 



(25) 



where we have defined a new scalelength £ including all the constants. Eq.(25) gives 
the potential for a pointlikc mass. In order to generalize this result to the case of 
an extended system, we may divide the system in infinitesimal mass elements dm 
and sum up the contributions to get the total potential. Assuming for the sake of 
simplicity spherical symmetry, the gravitational potential of an extended halo is 
thus : 



$(r) = - 



GM(r) 



1/3' 



4ttG / p(r r )r' 



1 1 



1/3' 



dr' 



(26) 



where M(r) and p{r) are respectively the halo mass and the density profile. The 
first term in Eq. (26) represents the contribution to the gravitational potential of the 
mass within the radius r, while the second one takes into account the mass outside 
this radius. The circular velocity depends on the force acting on the star orbiting 
at distance r from the halo centre and is thus only determined by the first term in 
Eq.(26) as a result of the Gauss theorem 47 . A straightforward generalization of the 
standard formula then gives in this case : 



v 2 c (r) = 



GM{r) 



i 2 ( r 



1/3 



(27) 



Comparing this result with the Newtonian formula v1(r) = GM(r)/r shows that 
the circular velocity is higher because of the additive term (r/^) 1 / 3 . Qualitatively, 
Eq.(27) shows that the fluid of strings effectively works as a fictitious source dis- 
tributed with a mass profile 2M(r)/3 x (r/^) 1 / 3 . Because of this additional effective 
source, it is possible to get the same value of v c (r) as in the Newtonian case with a 
smaller value of the CDM halo mass. Summarizing, introducing the fluid of strings 
modifies the gravitational potential in such a way that less CDM is necessary to get 
a given value of v c (r). 

The reduction of the quantity of CDM needed to fill the dark matter haloes is 
consistent with what is expected from our previous estimate of fib + flcDM in the 
ASDM model. Indicating with ft s CDM (^cdm) the CDM density parameter for the 
ASDM (ACDM) model, we may qualitatively write : 

''''CDM IV1 CDM 



CDM 



M N 



where M S CDM (Mg DM ) is the typical mass in CDM particles in a ASDM (ACDM) 
model. Naively speaking, in order to get the same value of the rotation curve at, for 
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instance, the virial radius f r vir , the halo mass in the ASDM case must be smaller 
than in the ACDM case by an amount that is of the order of magnitude of the 
correction term in Eq.(27) evaluated at r vir so that it is : 

^cdm ^ 2 / r V j r 

^CDM 3 V £ 

According to the results in Table 1, ^cdm/^cdm — 1/2 so that, from the above 
relation, we get £ ~ (4/3) 3 7\, ir . Such high values of £ also ensures that the gravita- 
tional potential is practically the same as the classical Newtonian one in the galactic 
regions dominated by the visible components where the rotation curve is well fitted 
by using the standard formulae. It is worth stressing, however, that this encouraging 
result needs to be further investigated by a careful fitting to the rotation curves of 
observed galaxies. This is outside the aim of this paper, but will be presented in a 
forthcoming work. 

5. Conclusions 

Shedding light on the dark side of the universe is a very difficult, but also very 
attractive challenge of modern cosmology. The nature and the fundamental proper- 
ties of the two main ingredients of the cosmic pie, namely the dark energy and the 
dark matter, are still substantially unknown and it is, indeed, this wide ignorance 
that justifies and motivates the impressive amount of theoretical models proposed 
to explain the observed astrophysical evidences. Moving in this framework, we have 
considered the dark matter as made out not only of massive dustlikc CDM particles, 
but also of a fluid of strings of cosmic origin with an equation of state w s = —1/3. 
Starting from this idea, we have considered four cosmological models comprising four 
components, namely dustlike baryons and CDM, fluid of strings and dark energy 
with constant barotropic factor w. Two of these four models (ASDM and QSDM) 
have a non vanishing fraction of dark matter in the form of a fluid of strings, while 
in two models (ACDM and ASDM) the energy budget is dominated by the A term. 
Our main results are briefly outlined as follows. 

(1) All the models are able to fit the data on the SNela Hubble diagram and the 
gas mass fraction in galaxy clusters with very good accuracy. In particular, it is 
remarkable that the total dark matter density parameter Qdm = ^cdm + 
is very well constrained and turns out to be the same in all models. When the 
assumption w = — 1 is relaxed, the dark energy barotropic factor is constrained 
to be in the region w < — 1 so that phantom like models are clearly preferred 
with a disturbing violation of the strong energy condition. It is worth noting 
that present data do not require phantom dark energy since they can be equally 



f The virial radius is defined such that the mean mass density within r v i r is S t h times the mean 
matter density p = S!m p cr it, with <5^ the critical overdensity for the gravitational collapse of 
density perturbations. 
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well fit by models with the cosmological constant A driving the accelerated 
expansion. Discriminating among the different possibilities will need a large 
sample of high redshift SNela such as those that should be available with the 
SNAP satellite mission 48 . 

(2) According to both the AIC and BIC, the QCDM model is statistically preferred 
over the other considered possibilities and this is not an unexpected result. 
However, this is obtained to the price of admitting phantom dark energy which 
is affected by serious theoretical difficulties. On the other hand, the ASDM 
model is preferred over the popular concordance ACDM scenario and is only 
slightly disfavoured with respect to the QCDM one. The good fit to the data 
and the graceful feature of avoiding to enter the realm of ghosts makes this 
model a good compromise between observations and theory and we therefore 
consider it as our final best choice. 

(3) The ASDM model predicts that a significant fraction (e ~ 59%) of the dark 
matter is made out by a fluid of strings so that the standard matter density 
parameter flf, + flc d m is only half of the fiducial value in the concordance sce- 
nario (0.15 vs 0.30). Since we may assume that the fluid of strings is massless 
(or nearly so), we should expect a corresponding decrease of the mass of galac- 
tic dark haloes. If the gravitational potential is still Newtonian, decreasing the 
CDM halo mass should lead to lower values of the circular velocity in the outer 
dark matter dominated regions of galaxies. This is not the case since, in the 
weak field limit, the ASDM model gives rise to a modification of the gravita- 
tional potential. As a result, the circular velocity due to a mass M(r) is higher 
than in the classical case so that less massive haloes are necessary to give the 
observed values of v c (r). Moreover, a very qualitative calculation suggests that 
the typical value of the scalelength over which deviations from Newtonian for- 
mulae cannot be neglected is sufficiently high that the inner luminous matter 
dominated rotation curve is unaltered. 

These encouraging results motivate further studies of the ASDM model. To this 
end, there are two different routes connected to two different features of the model 
which can be followed. 

First, because of its scaling with the redshift as p s oc (1 + z) 2 , that is intermedi- 
ate between that of CDM and that of A, a new era dominated by a fluid of strings 
is predicted in the expansion history of the universe. It is thus worth investigating 
how this imprints on the CMBR anisotropies in order to see whether the spectrum 
measured by WMAP is still accurately reproduced. To this regard, it is worth not- 
ing that the attempts recently made to constrain the cosmic strings contribution 
to the CMBR spectrum 49 may not be extended to our case since they refer to a 
network of cosmic strings rather than clouds of finite length strings. Less theoret- 
ically demanding, but more observationally ambitious is the possibility to test the 
proposed scenario on the basis of the transition redshift zt- As Table 1 shows, for 
the ASDM model, zt is significantly higher than in the other cases so that a model 
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independent estimate of this quantity could be a powerful discriminating tool. 

One of the most peculiar features of the ASDM model is the modified gravi- 
tational potential in Eq.(26) leading to the corrected circular velocity in Eq.(27). 
Having been obtained in the weak field limit, such correction should be tested at 
the scale of galaxies and clusters of galaxies thus offering the possibility to test the 
model at a very different level. To this aim, one should try fitting the rotation curve 
of spiral galaxies to see whether the problem of their flatness could be solved in this 
framework. Moreover, it is interesting to check how much the halo mass is reduced 
and to compare the reduction with respect to the classical Newtonian estimates 
with the decreasing of Vic dm obtained above. To this aim, low surface brightness 
(LSB) galaxies are ideal candidates since they arc likely dark matter dominated so 
that systematic uncertainties on the luminous matter modelling have only a minor 
impact on the fitting procedure. Moreover, the stellar mass - to - light ratio of LSB 
galaxies is well constrained so that we may fix this quantity thus decreasing the 
degeneracy among the other parameters. Useful samples of LSB galaxies with de- 
tailed measurements of the rotation curve are yet available (see, for instance, 50 ) 
so that this kind of test may be easily implemented. In this same framework, it is 
also interesting to consider the velocity dispersion curves in elliptical galaxies where 
recent studies seem to indicate a dark matter deficit 51 . 

Changing the gravitational potential does not only alter galaxies rotation curve, 
but also affects the clustering properties and thus leads to a different matter power 
spectrum. It is thus interesting to compare the predicted power spectrum with 
those measured from the SDSS galaxies in order to check the validity of the ASDM 
model. A similar comparison has been recently performed by Shirata et al. 52 for two 
phenomenological modifications of the law of gravity. We stress, however, that their 
approach is purely empirical and, furthermore, assumes that the universe can still 
be described at large scales with the ACDM model. Since in order to compute the 
power spectrum, one also needs the background Hubble parameter, it is important 
to use an expression for H(z) that is consistent with the proposed modification of 
gravity. For the ASDM model considered here, all the ingredients are at disposal 
so that a coherent calculation can be performed. It is worth noting that such a 
test is the only one capable of probing the model both at the galactic (through the 
gravitational potential) and cosmological (because of the use of H) scales at the 
same time. 

As a concluding general remark, we would like to stress the need for tackling the 
dark matter and dark energy problem together taking care of what is the effect at 
the galaxy scale of any modification of the fundamental properties of one of these 
two components. In our opionion, this could be a valid approach in elucidating the 
problems connected to the dark side of the universe. 
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